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I.  PLASMA  THEORY  AND  SIMULATION 


A.  Ion-ion  two-stream  mode  in  a  thermal-barrier  cell  ;  local  electrostatic  stability 
analysis 

V  V.  A.  Thomas  ( Dr.  W.  M.  Nevins )  '  '•'*  5 

\  a*  r,  , 

As  discussed  in  the  previous -QPRMII-  ,19ftt),  the  ion  distribution  function  in  the 
thermal-barrier  region  of  a  tandem  mirror  machine  oa^Be  unstable  to  electrostatic  and/or  elec¬ 
tromagnetic  two-stream  modes.  In  this  •QPR-we^examine^he  analytic  stability  in  an  infinite 
medium  of  those  electrostatic  ion-ion  two-stream  modes  which  have  thejr  wave  vector,  k, 
parallel  to  the  axial  magnetic  field  (this  is  equivalent  to  examining  1-d  'Idealized*'7 stability). 
The  point  of  our  doing  //(homogeneous  simulations  is  precisely  to  examine  any  non-infinite 
medium  effects  on  the  instability  in  a  thermal  barrier.  Any  unstable  ion-ion  two-stream  modes 
are  expected  to  have  wavelengths  of  orders  ;  since  the  thermal-barrier  cell  itself  is  much 
longer  than  this,  our  analytic  ("focal*)  approximation  will  provide  at  least  a  zeroth-order  esti¬ 
mate  of  instability  as  well  as  as  a  reference  from  which  to  measure  any  inhomogenity  effects 
recovered  later  by  our  particle  simulation  runs.  -i  \ 

Recently,  the  stability  of  colliding  Maxwellian  ion  beams  was  examined  by 
Foote  and  Kulsrud1.  Their  work  determined  the  stability  of  such  Maxwellian  ion  beams  to  both 
electrostatic  and  electromagnetic  modes,  assuming  an  infinite  homogeneous  plasma.  Because 
the  actual  thermal-barrier  cell  ion  beams  are  not  Maxwellian,  we  developed  a  code  to  calculate 
local  growth  rates  using  a  more  accurate  model  for  the  ion  distribution  function,  taken  from 
Cohen2. 

This  instability  requires,  roughly,  that  Tt/T,  >  3  and  that  the  relative  drift  between  the 
two  ion  beams  be  on  the  order  of  several  times  the  ion  thermal  velocity.  The  electron  contri¬ 
bution  to  the  analytic  dispersion  relation  is  modelled  in  our  work  as  just  the  real  quantity 
(tupt/kvtly2  -  (k\[))~2;  the  effect  of  this  electron  term  is  to  eliminate  the  ion-acoustic  insta¬ 
bility.  In  these  expressions,  it  is  important  to  understand  that  temperature,  T„  and  plasma  fre¬ 
quency,  Up,,  refer  just  to  the  thermal  electrons.  That  is,  if  a  "  hot "  electron  population  and  a  " 
thermal  ”  electron  population  are  both  present  at  the  same  location  with  Tyl0l»  Tltxemai,  then 
the  true  shielding  distance  is  determined  by  the  Debye  length  of  the  cooler  "  thermal  *  popula¬ 
tion.  This  is  an  important  consideration  because  if  the  thermal  electron  density  is  not  high 
enough,  then  the  ion  density  perturbations  are  not  shielded  sufficiently  to  prevent  instability. 
On  the  other  hand  a  high  density  of  hot  electrons  (which  are  magnetically  trapped)  is  necessary 
to  the  operation  of  the  thermal-barrier  cell:  these  maintain  the  large  dip  in  electrostatic  poten¬ 
tial  which  is  the  essence  of  the  thermal  barrier.  The  trade-off  between  these  two  competitive 
factors  is  be  explored  in  this  note. 

The  notation  of  Lontano  et  al. 3  will  be  followed  closely  here.  They  considered  essentially 
the  same  problem  except  that  their  distribution  function  did  not  include  thermal-barrier  trapped 
ions.  The  appropriate  dispersion  relation  for  our  problem  is 
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where  u  —  v„  is  the  ion  parallel  velocity,  and  where 

F(u)  - //,(▼)  5 (u-  k  v)  dv 

is  the  vn  'reduced*  ion  distribution  (i.e.,  /,  integrated  over  vj.  Here  k  is  a  unit  vector  in  the 
direction  of  the  wave  vector. 

A  standard  Nyquist  technique  is  used  to  determine  the  stability  of  this  system.  The  dielec¬ 
tric  function  Z(v,)  is  examined  as  v,  (representing  the  phase  velocity  <u/k)  is  varied  from  -<*> 
to  along  a  path  which  is  slightly  above  the  real  axis;  the  Nyquist  integration  path  is  then 


closed  by  a  half-circle  at  infinity  in  the  upper  half  u  plane.  The  number  of  unstable  modes  is 
given  by  the  number  of  times  the  function  Z(v,)  encircles  the  origin.  Figure  1,  taken  from 
Ref.  3,  shows  the  hodograph  for  this  integral.  In  our  case,  the  quantity  Im  (Z(v,»  is  seen  to 
equal  to  zero  at  only  five  locations  along  the  contour  path. 


R«u> 


FIG.  1.  Nyquist  hodograph  for  the  instability  analysis:  the  point  (1)  corresponds  to  v,  *  ±.  '0 
;  the  point  (2)  corresponds  to  »,  —  ±°°  and  (3)  corresponds  to  v,  —  0  .  The  solid  Im  (Z)  axis 
is  for  a  stable  distribution  and  the  dotted  Im  (Z)  axis  is  for  an  unstable  distribution.  Although 
the  figure  is  for  the  distribution  function  in  Ref.  3,  the  contour  topology  remains  the  same  for 
the  model  we  study,  taken  from  Ref.2  . 

These  five  points  are  v,  -  0  ,  v,  -  ±  «o  and  v,  —  ±  v0  where  these  v0  are  the  beam  velocities, 
where  the  reduced  ion  distribution  function  has  its  maxima.  From  the  figure  it  is  seen  that 
there  is  at  most  (me  unstable  mode  for  our  problem;  and  since  a *  —  a  +  /A  and  u  —  —  a  +  ib 
ate  both  solutions  to  our  dispersion  relation,  we  conclude  that  a  —  0  ,  and  therefore  that  the 
unstable  mode  is  purely  growing.  This  fact  facilitates  the  numerical  analysis  since  it  allows  us 
to  replace  «  by  />  in  our  numerical  search  for  the  unstable  roots.  A  simple  bisection  method 
can  then  be  used  to  determine  the  growth  rate  of  an  unstable  mode  for  a  given  test  problem. 

For  more  detailed  studies,  where  we  are  interested  in  more  than  just  the  unstable  solu¬ 
tions  to  the  dispersion  relation,  we  decided  to  use  a  packaged  complex  root  finder  utilizing 
Muller’s  method  (i.e.,  in  performing  these  calculations  we  do  not  substitute  iy  for  <■>  and  then 
use  bisection  to  locate  roots).  The  time  to  find  the  first  root  is  increased  but  all  subsequent 
solutions  have  good  initial  guesses  and  converge  rapidly.  This  also  allows  us  the  flexibility  to 
examine  all  solutions  to  the  dispersion  relation  in  cases  where  that  is  desired. 

We  have  checked  the  code  results  by  adjusting  the  parameters  of  our  model  distribution 
function  (taken  from  Ref.  2)  so  that  the  ion  beams  are  nearly  Maxwellian,  and  then  comparing 
the  stability  boundaries  recovered,  against  those  given  in  Ref.  1.  It  should  be  pointed  out, 
however,  that  the  distribution  function  for  the  actual  thermal-barrier  is  decidedly  non - 
Maxwellian,  especially  near  the  cell  midplane  where  u^/v/  v„  =  3  .  Note  also  that  the  work  of 
Ref.  1  neglects  the  quantity  (kXc,)2  compared  to  one:  although  this  is  valid  for  determining 
the  stability  boundary  (since  marginal  stability  occurs  for  kX^  equal  to  zero),  this  neglect  may 
not  always  be  valid  for  determining  the  growth  rate.  Stability  boundaries  using  our  ion  model 
could  well  differ  from  those  found  assuming  Maxwellian  ion  distribution  functions.  In  fact,  our 
initial  results  indicate  that  the  thermal-barrier  cell  using  our  model  ion  distribution  function  is 
significantly  more  stable  than  it  would  be  using  a  Maxwellian  ion  beam  model.  These  results 
will  be  presented  in  detail  in  future  publications. 
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B.  Alfven  Ion  Cyclotron  Instability  Particle  Simulation 

N.  Otani  (Dr.  W.  M.  Nevins  ) 

No  special  progress  to  report  this  quarter. 

C.  Plasma  Sheaths:  Electrostatic  l-d  Particle  Simulations 

S.  Rousset 

No  special  progress  to  report  this  quarter. 

The  title  will  change  to:  "Time- Dependent  Child- Langmuir 
Diode  Simulation”. 

D.  Plasma  Diode:  l-d  Vlasov  Simulation 

W.  Lawson  (  Dr.  T.  L.  Crystal ) 

No  special  progress  to  report  this  quarter. 

E.  Relativistic  Quantum  Conductivity  Tensor  for  a  Plasma 

Kwang-Youl  Kim  (  Dr.  T.  L.  Crystal ) 

No  special  progress  to  report  this  quarter;  further  work  will  be  given  in  a  complete  report. 
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n.  SECTION  n:  SUMMARY  OF  REPORTS,  TALKS,  PUBLICATIONS  AND 

VISITORS 

A.  ERL  REPORTS: 

(1)  Quasineutral  Hybrid  Simulation  of  Macroscopic  Plasma 

by  Douglas  S.  Hamed,  UCB/ERL  Mem.  M81/71  14  Sept  1981. 

(2)  Rotational  Instabilities  in  the  Field-Reversed  Theta  Pinch: 

Results  of  Hybrid  Simulations 

by  Douglas  S.  Hamed,  UCB/ERL  Mem.  M81/70  14  Sept  1981. 

(3)  Kink  Instabilities  in  Long  Ion  Layers 

by  Douglas  S.  Hamed,  UCB/ERL  Mem.  M81/73  16  Sept  1981. 

(4)  Stabilization  of  the  Lower  Hybrid  Drift  Instability  by  Resonant  Electrons 
by  Yu-Jiuan  Chen,  William  M.  Nevins  and  Charles  K.  Birdsall 
UCB/ERL  Mem.  M81/72  16  Sept  1981. 


B.  JOURNAL  ARTICLES: 

Strong  Ion  Ring  Equilibria  Formed  by  Iqjection  and  Intrinsic  Stochasticity  of  Orbits 
by  Alex  Friedman,  Jacques  Denavit,  and  R.  N.  Sudan 
Jour.  Comp.  Phys.  44  pp  104-133,  Nov  1981. 


University  of  California,  Berkeley 
Electrical  Engineering  &  Computer  Sciences 


Fall  1981  T.L. Crystal 

Fridays  3:00-5:00  PM  642-9672 

Room  299  Cory  Hall 

EECS  298-9:  PLASMA  THEORY  &  SIMULATION  SEMINAR 

Oct  2  Yu-Jiuan  Chen  and  Bill  Nevins,  LLNL 
Use  of  the  ZED  postprocessor. 

Oct  9  Tom  Crystal 

Kinetic  simulations  of  Tokamak  TEM 


Oct  23  Niels  Otani 

The  AIC  instability:  linearized  e-s  simulations. 


Oct  30  Mary  Hudson,  Cyndi  Cattell,  Ilan  Roth  &  Earl  Witt,  SSL 

Space  plasmas:  special  problems  being  studied  with  simulation 


Nov  6  Vince  Thomas 

Ballooning  instability  and  plasma  density  limitations  in 
the  Jovian  magnetosphere. 


Nov  13  AT  1:30  PM.  Gary  Smith,  LLNL 

Electron  microinstability  and  negative  reactors. 


Nov  20  Kwang-Youl  Kim 

Quiet-start  simulations  of  Rings  &  Spokes. 

Dec  4  Stephane  Rousset 

ID  diode  simulations:  progress  since  Langmuir's  solutions. 
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SUMMARY  0*  PROGRESS  FOR  FIRST  AND  SECOND  QUARTER,  1982 

This  is  a  summary  of  highlights  of  progress  made  by  our  group,  for  use  by  DOE  and  ONR 
contractors.  The  details  appear  in  our  Quarterly  Progress  Report. 

Our  group  uses  both  theory  and  simulation  as  tools  in  order  to  increase  the 
understanding  of  instabilities,  heating,  transport  and  other  phenomena  in  plasmas.  We 
also  work  on  the  improvement  of  simulation  both  theoretically  and  practically.  We  are 
the  undersigned  plus  five  graduate  students,  and  are  aided  by  professionals  from  the 
Lawrence  Livermore  National  Lab  and  the  Space  Sciences  Lab. 

A.  Ion-Ion  Two  stream  Mode  in  a  Thermal  Barrier  Cell 

The  model  uses  ion  particles  with  isothermal  Boltzmann  (T^cT^)  electrons  plus 
hot  ECRH  localized  electrons,  with  initial  values  obtained  from  a  Fokker-Planck  code. 

Ion  test  oarticles  advanced  in  the  equilibrium  fields  do  maintain  their  initial  distri¬ 
bution  function,  as  desired.  The  next  step  is  to  be  wholly  self  consistent. 

B.  Alfven  Ion  Cyclotron  Instability  Particle  Simulations  No  special  progress. 

C.  Tinie-Deoendent  Child-Langmuir  Diode  Simulation  The  particle  code  and  results  are 
in  an  ERL  reoort  to  be  issued  soon.  The  work  is  noteworthy,  with  especial  attention 
to  inlections,  apolied  potentials,  and  diagnostics,  for  one  species.  The  results  check 
Langmuir  verv  well. 

D.  Plasma  Diode: Id  Vlasov  Simulation  (GASBAG  Code)  This  electrostatic,  non-neutral. 

Id,  bounded  model  is  solved  by  integration  of  the  Vlasov  equation.  Jump  discontinuities 
in  f(x,v)  (e.g.,  due  to  accelerations  of  a  half-Maxwellian)  are  avoided  by  a  continuation 
procedure. 

E.  Weak  Monotonic  Double  Lavers  Analytic  solutions  for  two  DL's  (related  to  electron 
solitary  hole  and  to  ion  acoustic  solitary  hole)  are  presented  in  an  ERL  report 

to  be  Issued  soon. 

F.  Plasma  Sheath  Formation.  Ion  Acceleration  and  Fluctuations  in  Steady  State  At  a 
floating  orobe  large  Up  oscillations  in  potential  are  observed  (in  simulation) ,  with 
ion  acceleration  to  o(vsoun(j  (Ko<<2)  at  the  wall.  (Abstract) 

G.  Sheath  Formation  and  Fluctuations  with  Dynamic  Electrons  and  Ions  Similar  to  F, 
with  results  for  1  <  mi/me  <1600  for  Te»  T^  and  for  0.1  <  Tj/Te  <  64  for  mi/me*400, 
showing  the  transient  formation  of  a  sheath. 

H.  Electron  Diode  Dynamics;  Limiting  Currents;  Plasma  Diodes  A  double  layer  is  viewed 
as  a  virtual  cathode  next  to  a  virtual  anode,  with  expectations  of  both  current  (flux) 
limiting  and  oscillations. 

Code  Development:  The  postprocessor  ZED  is  now  available  on  the  MFE-CRAY. 

One  ERL  report  was  issued  in  this  period  (already  distributed,  Birdsall  sheath  survev) : 
three  talks  were  presented  at  meetings;  visits  were  made  to  four  foreign  labs. 
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SUMMARY  OF  PROGRESS  FOR  FOURTH  QUARTER 
1981 

This  is  a  suomary  of  highlights  of  progress  made  by  our  group,  for  use  by  DOE  and  ONR 
contractors.  The  details  appear  in  our  Quarterly  Progress  Report. 

Our  group  uses  both  theory  and  simulation  as  tools  in  order  to  Increase  the 
understanding  of  instabilities,  heating,  transport  and  other  phenomena  in  plasmas. 

We  also  work  on  the  improvement  of  simulation  both  theoretically  and  practicallv. 

Our  group  consists  of  Prof.  Birdsall,  Dr.  Crystal  and  five  graduate  students. 

We  are  aided  by  professionals  from  the  Lawrence  Livermore  National  Lab  and  the 
Space  Sciences  Lab. 

A.  Ion-Ion  Two-Stream  Mode  in  a  Thermal  Barrier  Cell;  Local  Electrostatic 
Stability  Analysis.  The  system  studied  has  opposing  ion  streams,  with  fj(v) 

obtained  from  Fokker-Planck  thermal  barrier  studies  of  R.  H.  Cohen  (1981), 
and  warm  electrons.  The  Nyquist  technique  was  used  to  locate  regions  of  instability 
and  a  root  solver  used  for  detailed  roots.  Initial  results  show  that  the  thermal 
barrier  cell  is  significantly  more  stable  than  it  would  be  with  a  Maxwellian  ion 
beam  model. 

Research  on  the  following  topics  continue,  with  progress  to  be  reported  later: 

B.  Alfven  Ion  Cyclotron  Instability  Particle  Simulation 

C.  Plasma  Sheaths:  Electrostatic  lu  Particle  Simulations 

D.  Plasma  Diode:  Id  Vlasov  Simulation 

E.  Relativistic  Quantum  Conductivity  Tensor  for  a  Plasma 


Four  ERL  reports  were  Issued  and  distributed  this  quarter  and  one  journal 
article  published.  Two  reprints  are  distributed  with  this  Q*»R;  both  by 
Douglas  S.  Harned: 

"Quasineutral  Hybrid  Simulation  of  Macroscopic  Plasma  Phenomena",  Jour.  Comput  Phvs. 
Phys.  47,  pp  452-462  Sept.  1982. 

"Kink  Instabilities  in  Long  Ion  Layers"  Phys.  FI.  25  pp  1915-1921,  Oct.  1982 
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